Spin Configurations and Activation Gaps of the Quantum Hall States in Graphene 
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We report on our accurate evaluation of spin polarizations of the ground state and particle-hole 
gaps for partially-filled lowest Landau level, observed in recent experiments on graphene subjected 
to ultra-high magnetic fields. We find that inter-electron interactions are important at these filling 
factors, characterized by the non-degenerate ground states and large particle-hole gaps at infinitely 
large wave vectors. The gaps are the largest for the quantum Hall states in the second Landau level. 
The weak appearence of quantum Hall plateaus in the experiments for certain filling factors in the 
second Landau level however indicates that at these filling factors, the system has a soft mode at a 
finite wave vector due presumably to the presence of disorder. 
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In recent experimental work, an atomically thin, two- 
dimensional (2D) sheet of graphite, called graphene has 
been extracted by micromechanical cleavage A two- 
dimensional electron system in graphene exhibits many 
remarkable properties. In the band structure calculations 
where electrons are treated as hopping on a hexagonal 
lattice [H, one finds a unique linear (relativistic type) 
energy dispersion near the corners of the first Brillouin 
zone where the conduction and valence bands meet. As a 
consequence, the low-energy excitations follow the Dirac- 
Weyl equations for massless relativistic particles In 
an external magnetic field, the spectrum develops into 
Landau levels, each of which approximately fourfold de- 
generate [3, [f|. Recent discovery of the quantum Hall 
effect in graphene 0, 0, @l nas resulted in intense ac- 
tivities (JllO; El to unravel the electronic properties of 
graphene that are distinctly different from the conven- 
tional (or, as commonly called the 'non-relativistic') 2D 
electron systems in semiconductor structures. 

Graphene has a honeycomb lattice structure of sp 2 car- 
bon atoms, with two atoms A and B per unit cell (i.e., 
a two-dimensional triangular Bravais lattice with a basis 
of two atoms). Each atom is tied with its three nearest 
neighbors via strong a bonds that are in the same plane 
with angles of 120°. The it orbit (2p z ) of each atom is per- 
pendicular to the plane and overlaps with the n orbitals 
of the neighboring atoms that results in the delocalized 
7r and 7r* bands. There is only one electron in each 7r 
orbit and the Fermi energy is located between the tt and 
7r* bands [T^]. The dynamics of electrons in graphene is 
described by a nearest-neighbor tight-binding model Q 
that describes the hopping of electrons between the 2p z 
carbon orbitals. The first Brillouin zone is hexagonal and 
at two of its incquivalent corners (the K and K' points) 
the conduction and valence bands meet. Graphene is of- 
ten described as a two valley (K and K') zero-gap semi- 
conductor. Near these two points (the so-called Dirac 
points), the electrons have a relativistic-like dispersion 
relation, e k = ±hv\k\ and obey the Dirac- Weyl equations 



for massless fermions. At the vanishing gate voltage, the 
system is half-filled and the Fermi level lies at the Dirac 
points. 

In graphene, discovery of the Hall conductivity a xy = 
v G e 2 /h = 2(2n+l)e 2 //i, where n is an integer has indeed 
confirmed the Dirac nature of the electrons in that sys- 
tem [f| . This result has been explained as an outcome of 
four times the quantum Hall conductivity (n + \)e 2 /h for 
each Dirac fcrmion, that is consistent with the four-fold 
degeneracy (two spins and two valleys.or flavors) of the 
Landau levels for the Dirac fermions |8[. In this paper, 
we however focus our attention to the extreme quantum 
limit. In the very high magnetic field regime (up to 45 
Tesla), Zhang et al. [7[ discovered a new set of quantum 
Hall states at filling factors v G = 0,±I,±4. Observa- 
tion of these states clearly indicated that the four-fold 
degeneracy of the observed states at v Q = ±4(n + |) 
has been lifted in the presence of these ultra-high fields. 
According to Nomura and MacDonald (HI, appearence 
of the half- and quarter-filled lowest Landau levels are a 
consequence of quantum Hall ferromagnetism induced by 
the electron-electron interactions. The absence of certain 
Hall plateaus, such as at v G = ±3, ±5 are then attributed 
to the presence of disorder in the system. 



Alicea and Fisher [I4| performed a comprehensive ex- 
ploration of alternative mechanisms for these intermedi- 
ate quantum Hall states in clean and dirty graphene sam- 
ples. In particular, they looked for mechanisms which ex- 
plicitly could break the inherent symmetries of graphene. 
They started from the lattice scale with a Hubbard type 
model introducing inter- and intra-sublattice interactions 
to break the underlying SU(4) symmetry of the ideal 
graphene lattice. In the ferromagnetic (clean samples) 
regime they found the system to favor flavor polariza- 
tions at filling factors v G = ±1. The strength of interac- 
tions relative to the Zceman coupling determines whether 
the system is polarized or unpolarizcd at v G = 0. In 
the higher Landau levels, these authors predicted that 



2 



skyrmions would provide the minimum energy charge ex- 
citations. At the filling factor v G = ±4, interactions lead 
to a spin-polarized ground state whereas the states at 
Vq = ±3, ±5 break spontaneously the flavor symmetry. 
In the paramagnetic (dirty samples) regime the quantum 
Hall states may be possible at v G — 0, ±1, ±4 but not at 

Vq = ±3, ±5. 

The question of the absence of QH steps at v G = 
±3, ±5 in the experiment has been raised by Fuchs and 
Lederer [l5|, who, as opposed to the models described 
above, suggest that electron-electron interactions play 
no role in the formation of integer quantum Hall states. 
Instead they propose that the inversion symmetry of 
graphene is spontancosly broken by a magnetic field de- 
pendent Pierls distortion. In the lowest Landau level this 
leads to the splitting of the valley degeneracy. This com- 
bined with the Zeeman splitting explains the existense of 
the observed Hall plateaus at intermediate filling factors 
Vq = 0, ±1. In the higher Landau levels, however, the 
lattice distortion cannot lift the valley degeneracy. This 
leaves only the Zeeman splitting leading to intermediate 



plateaus only at v ( 
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We however believe that inter-electron interactions 
play an important role for the QH states at partially- 
filled Vq = 1 that do or do not appear as QH steps in 
the experiments. In this paper, we set out to accurately 
determine the spin configurations and activation gaps of 
various integer filling factors that we expect to shed light 
on the nature of these states as well as provide a better 
understanding of the source of lifting of degeneracy for 
these states. In conventional quantum Hall states at par- 
tially filled Landau levels the electron-electron Coulomb 
interaction opens up an excitation gap between the dif- 
ferent, otherwise degenerate eigenstates of the total spin. 
In graphene, an additional mechanism is necessary to lift 
also the degeneracy between the eigenstates of the total 
flavor pscudospin. Motivated by the lattice scale studies 
by Alicea and Fisher [lij we introduce into our Hamil- 
tonian a term which on the lowest Landau level maps to 
a charge imbalance between the sublattices. We found 
that the ground states of v G = ±1,±2,±3, etc. are 
non-degenerate, and they have large particle-hole gaps 
at asymptotically large wave vectors - the quasiparticle- 
quasihole gap for the fractional filling factors 1(1, 
which is in fact, the activation gap. Interestingly, that 
gap is much larger for the QH states at the n = 1 Landau 
level, as compared to that in other Landau levels, which 
is in line with earlier analysis of the effects of Coulomb 
in graphene. However, in the experi- 
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interactions ^ 

ments Q reported as yet, the absence (or weak presence) 
of the QH states at v G = ±3, ±5, etc. would imply that 
the system has a soft mode at a finite wave vector, which 
might arise due to disorder [l4[ in the system. 

In our model the electrons reside on the surface of 



a torus formed from a rectangular cell of sides L x and 
L y by imposing periodic boundary conditions in both 
x and y directions. Choosing the Landau gauge B = 
V x (0, Bx, 0) T , we get the two-component wavefunctions 
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for the valley K and 
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-isgn(A) 4> n -l,m{z,y) 

4>nm{x, y) 



4>nm(x,y) 
-isgn(A)0 n _i, m (x,y) 



for the valley K' . Here n is the Landau level index and A 
is its energy ±\/2n in units of j/i, where I is the magnetic 
length defined as £ 2 = ch/eB. In the case of the lowest 
Landau level, n = 0, the components with n — 1 in the 
above functions are taken to be zero. The component 
wavefunctions (f> nm are the Laundau wavefunctions for 
the non-relativistic systems, 
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x H n ((x + qL x -X m )/£) 

where the centers X m of the harmonic oscillations are 
given by 
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Due to the periodic boundary conditions, the quantum 
number m yielding non-equivalent wavefunctions is re- 
stricted to the values m = 0, 1, . . . , M — 1, where M is 
the degeneracy of the Landau level and also equal to the 
number of flux quanta through the rectangle, i.e.. 
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We diagonalizc the interaction Hamiltonian 
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1^3 = 1 



in the basis 
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0,1, 
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of the four-component vectors. When appended with the 
spin degree of freedom we have effectively a set of eight- 
component vectors making the base prohibitively large 
for practical computations. We can, however decrease 
the number of basis states by noting first that there is no 
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spin dependent term in our Hamiltonian (Zeeman ener- 
gies can always be added afterwards) so the total S z is a 
good quantum number. If we fix S z to its minimum value 
(0 or 1/2) the base is significantly reduced but it still al- 
lows us to recover all possible eigenstates of the total S 2 . 
Secondly because the Coulomb interaction preserves the 
quantity m = Wj mod M, performing diagonalization 
separately for each m reduces the size of the basis further 
approximately to the 1/M-th part. 

We are still left with the flavor symmetry associated 
with the valleys K and K' leading to the degenerate 
ground states. Because on the lowest Landau level in 
the valley K (K 1 ) only the sublattice B (A) component 
of the wavefunction differs from zero the electron occu- 
pation of the valleys maps directly to the occupation of 
the sublattices. We also know that the electrostatic en- 
ergy is minimized if the electrons reside only in one of the 
sublattices (see [HI)- We utilize this fact by introducing 
the symmetry breaking term 



plicit expression 
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where K is the pseudospin operator which is identical 
to the spin operators S but works in the flavor (or val- 
ley) space. Introduction of this term in the Hamiltonian 
would imply that either one distributes the charge evenly 
(v KK > 0, Ki and Kj antiparallel) between both A and 
B sublattices or favor one of the sublattices (v KK < 0, 
Ki and Kj parallel) over the other. The latter option ap- 
plies then only to the lowest Landau level. In the higher 
Landau levels this interpretation is not valid since the 
wavefunctions there have equal weight on both sublat- 
tices and the pseudospin operator cannot redistribute the 
charges. Therefore, the purpose of the K ■ K term in the 
present case is only to lift the valley degeneracy. We 
may also speculate (consistently with the experiments) 
that this kind of symmetry-breaking term exists only in 
the lowest Landau level while in the higher levels a new 
(yet unknown) mechanism may operate with increasing 
magnetic fields. 

To construct the interaction Hamiltonian we need to 
evaluate the two-body matrix elements 
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where J\f is the normalization constant and we have de- 
fined the quantities Z as 



2AM 



(ia + sgn(n — n')(3X) 



n-n'| / n 1^2 , a 2\2 s 



Dim(ra,n') V. AM 
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A straightforward but lengthy algebra leads to the ex- —0.001 in units of Coulomb energy for v 



In the summation (with prime) the term with a = (3 = 
is excluded and the index (3 is restricted to the values 
j3 = kM + ni4 — mi; k = 0, ±1, ±2, . . .. For the nonzero 
matrix elements the quantum numbers m; must satisfy 
the conservation law mi + m 2 = ms + mod M. 

We have evaluated the total spin S, the ground state 
energy e coul (the Landau level energy is not included) 
and the activation gap (defined as thepositive disconti- 
nuity of the chemical potential fj.) [lil [20l | , a measurable 
quantity 

Ag = fj,+ - H- 

w -2A^ e e(^ G ) + N e e{v+) + N e e{v G ) 

for Landau levels n = 0,1,2,3 and for filling factors 
v L = N e /M with N e electrons and M flux quanta in 
the Landau level at issue. Here +(— ) refers to the case 
when a particle is added to (subtracted from) the system, 
and £{vq) is the total energy of the system at the filling 
factor v G . The filling factor v G for graphene is related to 
the conventional (i.e., non-relativistic) filling factor [l6| 
v L + An as 

Vq = Vj^ + 4n — 2. 

Our results for various integer filling factors are sum- 
marized in Table I and Table II for systems containing 
N e = 4 — 12 electrons. Clearly, the total spin is S = N e /2 
for the filling factor v L = 1 and 5 = for the filling factor 

v L = 2 that are independent of the Landau levels. The 
ground state (correlation) energies and the activation 
gaps are presented in Table I for a system with four flux 
quanta and in Table II in the case of six flux quanta. Due 
to the particle- hole symmetry the states with v G = 4n±l 
arc equivalent and the ones with v G = An + 2 correspond 
to completely filled Landau levels. Therefore they are not 
shown here. In our calculations, we have used the value 

It should 
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TABLE I: The Coulomb energy, the activation gap, and the 
ground state spin for four and eight electrons and four flux 
quanta. Here n is the Landau level index, e CO ul, and A are in 
units of the Coulomb energy (e 2 /d, e — 1). 
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0.4720 






TABLE II: Same as in Table I, but for six and twelve electrons 
and six flux quanta. 
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0.9810 


3 


2 


8 
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11 
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12 


-0.5656 


0.5004 






be pointed out, however, that the results presented in the 
tables are not sensitive to the magnitude of this param- 
eter: in fact, an infinitesimal value would be sufficient to 
yield, e.g., the total spins. 

The non-degenerate ground states (at zero wave vec- 
tor) for these filling factors, as well as large activation gap 
that appears at the very large wave vector (lfl . 17 1 implies 
that the system is incompressible at these filling factors. 
Interestingly, as the results indicate, the activation gap 
is always the largest for n = 1 Landau level. As men- 
tioned above, this is in line with earlier works [3, [HI, 
where electron correlations were found to be dominant 
for the n = 1 Landau level than for all the other Lan- 
dau levels, a significant shift from the conventional 2D 
electron systems. In the present case that would mean 
that the quantum Hall states should be most stable for 
v G = ±3, ±5, ... but in fact, they are conspicuously ab- 



sent in the reported experimental results [7j. Given the 
incompressible nature of these states for the large wave 
vector and for the zero wave vector as found above, we 
can only conclude from this apparent contradiction that 
these states must exhibit soft modes at a finite wave vec- 



tor, perhaps due to the presence of disorder ljl bi, or 
maybe due to a mechanism of the type proposed in [15j 
that inhibits the opening of the gap at small wave vec- 
tors. 

We would like to thank Jason Alicea for helpful dis- 
cussions. The work has been supported by the Canada 
Research Chair Program, a Canadian Foundation for In- 
novation Grant and the NSERC Discovery grant. 



[i] 
[i] 



[9 

[10 

[11 
[12 



[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 



95, 
96, 

183 



Electronic address: tapash@physics.umanitoba.ca 
K.S. Novoselov, et al., PNAS 102, 10451 (2005); Science 
306, 666 (2004); Y. Zhang, et al., Phys. Rev. Lett. 94, 
176803 (2005). 

P.R. Wallace, Phys. Rev. 71, 622 (1947). 

T. Ando, in Nano-Physics & Bio- Electronics: A New 

Odyssey, edited by T. Chakraborty, F. Peeters, and U. 

Sivan (Elsevier, Amsterdam, 2002), Chap. 1. 

J.W. McClure, Phys. Rev. 104, 666 (1956); R.R. Haering 

and P.R. Wallace, J. Phys. Chem. Solids 3, 253 (1957). 

Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 

K.S. Novoselov, et al., Nature 438, 197 (2005); Y. Zhang, 

Y.-W. Tan, H.L. Stormer, and P. Kim, ibid. 438, 201 

(2005) . 

Y. Zhang, et al., Phys. Rev. Lett. 96, 136806 (2006). 
V.P. Gusynin and S.G. Sharapov, Phys. Rev. Lett. 
146801 (2005); E. McCann and V.I. Fal'ko, ibid. 
086805 (2006). 

A.K. Geim and K.S. Novoselov, Nature Mat. 6, 
(2007). 

M.I. Katsnelson, Materials Today 10, 20 (2007). 
T. Chakraborty, Physics in Canada 62, 351 (2006). 
R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical 
Properties of Carbon nanotubes (Imperial College Press, 
London, 1998). 

K. Nomura and A.H. MacDonald, Phys. Rev. Lett. 96, 
256602 (2006). 

J. Alicea and M.P.A. Fisher, Phys. Rev. B 74, 075422 

(2006) . 

J.-N. Fuchs and P. Lederer, Phys. Rev. Lett. 98, 016803 

(2007) . 

T. Chakraborty and P. Pietilainen, The Quantum Hall 

Effects (Springer, Heidelberg, 1995), 2nd edition. 

C. Kallin and B.I. Halperin, Phys. Rev. B 30, 5655 

(1984). 

M.O. Goerbig, R. Moessner, and B. Doucot, Phys. Rev. 
B 74, 161407 (2006). 

V.M. Apalkov and T. Chakraborty, Phys. Rev. Lett. 97, 
126801 (2006). 

T. Chakraborty, P. Pietilainen, and F.C. Zhang, Phys. 
Rev. Lett. 57, 130 (1986). 



